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Abstract 

We consider formation of composite strings and domain walls as a 
result of fusion of two elementary objects (elementary strings in the 
first case and elementary walls in the second) located at a distance 
from each other. The tension of the composite object T2 is assumed 
to be less than twice the tension of the elementary object Xi, so that 
bound states are possible. If in the initial state the distance d between 
the fusing strings or walls is much larger than their thickness and 
satisfies the conditions T\d 2 3> 1 (in the string case) and T\d? 3> 1 (in 
the wall case), the problem can be fully solved quasiclassically. The 
fusion probability is determined by the first, "under the barrier" stage 
of the process. We find the bounce configuration and its extremal 
action Sb- In the wall problem e~ Se gives the fusion probability 
per unit time per unit area. In the string case, due to a logarithmic 
infrared divergence, the problem is well formulated only for finite- 
length strings. The fusion probability per unit time can be found in 
the limit in which the string length is much larger than the distance 
between two merging strings. 



1 Introduction 



In this paper we will consider two problems of practical interest which arise in 
various settings, and can be solved purely quasiclassically. The formulation of 
the problems, and their solution, is very general. They refer to restructuring 
of solitonic objects (or branes) supported in various field theories. We were 
motivated by a specific problem that arose in [1], but here we will give a 
general discussion, and find a generic solution, so that our results can be 
used in all similar situations. 

The first problem is about strings (flux tubes). Suppose we have two types 
of strings: "elementary" strings with tension Ti, and a composite string with 
tension T 2 . We assume that the composite string is a bound state, i.e. 

T 2 - 2Ti < . (1) 

By composite we mean that there is a conserved "charge" Q, and Q — 1 for 
the elementary string while Q = 2 for the composite one. The composite 
string can form as a result of a fusion of two elementary ones. 

One can consider two parallel strings at a distance d from each other. 
The parameter d is assumed to be much larger than the string thickness. 
Quantum fluctuations of strings can result in a configuration with two ele- 
mentary strings forming a composite one in the middle (see Fig. [1]). A crucial 
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Figure 1: Two elementary strings merge into a composite one. Once the size £ of the 
merged segment reaches (in a quantum tunneling process) its critical value £ c , further 
expansion of its size becomes classical. 

characteristic of the process is a critical size l c of the merged segment. It is 
determined by the balance of two energies: the one gained due to the fact 
that T 2 — 2Ti < (the energy gain is t c (2Ti — T 2 )), and the one lost due 
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to elongation of the stringSQ The first stage of dynamics producing a size 
£ c merged segment occurs as a quantum tunnelingjj which can be described 
in the Euclidean space-time. Once the critical-size segment is attained, its 
further expansion proceeds as a purely classical process, with positive energy 
release and an accelerating expansion of the merged segment. The fusion 
probability is determined by the quantum tunneling stage. 
The action corresponding to such fusion is large, provided 



As was mentioned, we gain energy in the central domain because T2 < 2T%. 
We loose outside the central domain because of the string bending needed 
to match the asymptotic boundary conditions. An extremal (in fact, maxi- 
mal) value of the action must exist. It is realized on a classical solution in 
Euclidean. This is a maximum with respect to the size variation i. This is 
the only instability, and is the usual one that gives the vacuum decay. 

The problem is similar, in a sense, to that of metastable vacuum decay [2, 
3] (for reviews see [4,5]), but with an important difference (see Appendix A). 
In the false vacuum decay the energy balance is achieved between a bubble of 
a genuine vacuum (gain) versus the potential energy of its surface (loss). The 
Euclidean solution is provided by a bounce configuration. In the problem at 
hand, in which the three-string junction is assumed to carry no energy, the 
barrier is not due to the potential energy, but is rather associated with the 
kinetic energy term in the string Lagrangian. However, this is just a technical 
difference. A critical field configuration extremizing the Euclidean action still 
exists, and we will find it in the limit when the string thickness is negligible 
compared to the interstring distance d, see Fig. [TJ This is an analog of the 
standard bounce [41 n 



In the formulation of the string fusion problem there are infrared sub- 
tleties related to the tails of the strings. We can define the problem in a 

^n the present work it is assumed that the three-string junction mass is negligible. 

2 This is due to the fact that at £ < l c the energy gain is less than the energy loss; 
therefore the system under consideration tunnels under a barrier. 

3 We will apply this term, bounce, to the extremal Euclidean string and wall configu- 
rations in the problems to be discussed below. 



Ti d 2 > 1 . 



(2) 
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finite box, of length L, and compute, with exponential accuracy, the prob- 
ability per unit of time of this fusion process T(L). This is the content of 
Section [2j 

The result of the above computation immensely simplifies if one calculates 
the exponent in a logarithmic approximation. In this approximation it turns 
out possible (in the limit 27\ — T 2 <C Ti) to generalize the analysis of the 
parallel string fusion to cover the case of nonparallel strings. This problem 
will be addressed in Sect. [3j 

The third problem we will deal with, in Section HI is similar in nature, but 
it refers to parallel domain walls, rather than strings. Adding an extra di- 
mension to the solitonic objects to be fused has a crucial effect. The infrared 
problem we had to deal with in the case of the string fusion now disappears, 
even in the infinite volume. Then, we can readily calculate the fusion prob- 
ability per unit time and per unit area, with the exponential accuracy. This 
is done in Sect. HI In Sect. Owe deal with the wall merger at strong binding. 
An instructive example of the domain wall fusion in super- Yang-Mills theory 
in considered in Sect. [61 

Is not difficult to generalize our analysis to branes with arbitrary p spatial 
dimensions, usually called p-branes, in a space-time with D + 1 dimensions. 
We outline this, and summarize our results in Sect. CD 

In summary, our solution of the string/wall fusion problem is general and 
independent of dynamics of the underlying microscopic theory provided the 
following assumptions are met: (i) 5/d — > where 5 is the string or wall 
thickness; (ii) Tid 2 (in the string case) or Tid 3 (in the wall case) 1; (iii) 
the three-string (or three- wall) junction contribution to the extremal action is 
negligible. The latter condition is met in many instances of practical interest. 
In the string case we must also assume that L ^> d. At weak binding the 
constraint on d softens; it is sufficient to require Tid 3 y / ^^r 3> 1. 

2 Parallel Strings 

To compute the decay probability, it is convenient to Wick-rotate the time 
direction. Then in the Euclidean space-time we have a problem of two static 
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2-branes, which can fuse due to quantum fluctuations]^! The Euclidean string 
action is the string tension multiplied by the area of the branes. 

We want to find a bounce solution, which corresponds to two surfaces 
at asymptotic distance d, and an interior "bubble" in which they overlap to 
form a bound state (Fig. [2]). The tunneling rate is then determined by the 
difference Sb between the (Euclidean) action on the bounce configuration 
and that on the trivial configuration, with two flat world sheets for each 
string, one at z = d/2 and another at z = —d/2, 

r = C exp(-5s) , (3) 

where C is a pre-exponential factor. Getting this factor requires calculation 

I 1 st string worldslieet 




Figure 2: World surface for two elementary strings forming a composite one. The 
Euclidean time is denoted by r. 

of the path integrals over fluctuations around the bounce solution, as well 
as around the trivial flat world-sheet configurations. This issue will not be 
addressed here. In what follows we will discuss only the exponential factor 
determined by classical solutions to the Euclidean equations of motion. 

One can see, however, that for infinite strings a bounce solution does 
not exist. This is due to an infrared peculiarity of two-dimensional surfaces. 
Whenever we pull such a surface in the perpendicular direction, it never 
becomes flat asymptotically. Its asymptotic behavior, from the solution to 

In Euclidean space one can view these fluctuations as thermal. 
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the Laplace equation, is always logarithmic. Thus, we have no chance to 
recover the required boundary condition, that at x, r — > oo the z coordinate 
of the surface tends to ±d/2. 

This infrared behavior can be regularized provided we assume that z = 
±d/2 is achieved at some finite distance in the {x, r} plane. The most 
physically transparent regularization of the Euclidean version of the problem 
consists of a strip (infinite in the r direction), with the boundary conditions 
z = ± d/2 implemented at its edges (Fig. [3]). We parametrize the coordinates 
as x, t and z, where x = =F-L/2 present two edges of the strip, r corresponds 
to the Euclidean time, and z = ±d/2 are the vertical locations of the two 
parallel strings, so that the boundary condition for the bounce configuration 
sought for is as follows: at x = ±L/2 the value of z is fixed at +d/2 for one 
string and at — d/2 for the other. We use x and r to parametrize the brane, 
and z = f(x,r) determines the height of the branes. In fact, it is sufficient 
to consider only the upper side of the picture since it is symmetric under 
reflection z — > —z. This will be referred to as a "strip" boundary condition. 

Below we will find that Sb depends on L only through In L. Aiming at log- 
arithmic accuracy (i.e. keeping InL and omitting nonlogarithmic constants 
assuming InL to be large), we can replace the strip boundary conditions by 
much simpler ones, to be referred to as "round" boundary conditions (Fig. [3]). 
The round boundary condition is convenient for two reasons. First, it will 
help us to calibrate our solution. Second, the results obtained with the round 
boundary condition are useful in extending the problem to the case of non- 
parallel strings. The problem with the round boundary conditions is that, by 
itself, it has no Minkowski physical interpretation. In the Euclidean space, 
instead, it is just a problem of fusion due to thermal fluctuations, with the 
position of the 2-branes fixed at the circle. 

2.1 String fusion with round boundary conditions 

In this section we will require z(R) = ±d/2 where 

R > d (4) 

is assumed (for the definition of R see Fig. [3]). In the leading logarithmic 
(in R) approximation the problem of merging of two parallel strings can be 
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Figure 3: Boundary conditions for the fusing string worldsheet. Solid lines: "strip" 
boundary conditions. Dashed line: "round" boundary condition. 

solved for arbitrary relation between the tensions T 2 and T\ as long as the 
merger is possible, i.e. T 2 < 2T\. We do not have to require 27\ - T 2 <C 
Thus, in this section we lift this constraint. 

Thus, we replaced the strip space-time boundary for the world sheets by 
a disk of a large radius R ~ L/2, so that a bounce centered at the origin (i.e. 
x = and r = 0) is 0(2) axially symmetric and is described by a function 
z(r) where 

r = V x 2 + t 2 . (5) 

The slice of the solution z(r) passing through the x = r = line is shown 
in Fig. [H (where t c = 2r c and L must be replaced by 2R). The contribution 
of each string's world sheet to the action for such a centrally symmetric 
configuration is given by the integral 

S = 2nT J rdr y/l + z' 2 (6) 

with an appropriate tension T. Here z' = dz/dr, and the integrand represents 
the area of the circular element of the surface. 
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The central part of the bounce configuration, a disk of radius r c located 
at z = 0, is filled by the string with tension T 2 . The profile z(r) for each 
of the strings with tension T\ is determined by the equations of motion, 
which extremize the surface area of two world sheets. The difference Sb of 
the action on the bounce and that on the trivial configuration can be thus 
written as 

S B = Tv(T 2 -2T 1 )r 2 c + A7rT 1 f rdr (y/l + z' 2 - l) , (7) 

where z(r) stands for the vertical profile of one of the two world sheets (for 
definiteness we consider the upper one) and the contribution of the other 
simply doubles the coefficient in front of the integral in Eq. (JTj). 

We find that the simplest way to analyze the solution for z(r) is using an 
"integral of motion" , which follows from the symmetry under z translation, 
which we call r^. 

r z' 

r ■ (8) 



The left-hand side is independent on z, which, in fact, tells us that the 
vertical component of the capillarity force acting on any horizontal section 
of the film is constant. The relation between the constant and the radius 
r c of the bounce is found from the condition of equilibrium of the boundary 
of the disk, where the string T 2 bifurcates into two strings 7\. This condition 
is that the net horizontal force at the boundary vanishes, 



2Ti 



,12 



T2 ■ (9) 



After setting r = r c in Eq. ([8]) and eliminating z' \ r =r c from Eqs. ([8]) and 
one readily finds the following relation: 



ro = M/!-J|- (10) 
The solution to the equation of motion ([8]) satisfying the boundary condition 
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z(R) = d/2 at R > r c has the form§ 



ro In ■ 



2R 



d 

+ 2' 



The parameter r c can be determined from the condition z(r c ) = in terms 
of d and R. To this end we substitute Eq. (jTOl) in 



In 



2R 



d 



rg 2r 



(12) 



which can be solved numerically. Figure H] presents r c /d as a function of R/d 
at a representative value of T 2 /(2Ti) = 0.95. (A matching with Eq. (fl9~I) 
below starts emerging at the right edge of the plot.) 
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Figure 4: The plot of r c /d vs. R/d at ^ = 0.95 (i.e. <5 » 0.31, see Eq. ([T7| ). 

Given the solution (fTTj) it is not difficult to find the action ([7]). It turns 
out that the action Sb has the simplest form being expressed in terms of Tq 

D To be more exact, in Eq. (JTUJ) z(R) — d/2 up to terms O . As we will see 

shortly, roughly speaking, r$ ~ dj [2 ln(i?/d)]. Hence, up to logarithms, the relative error 
is O ( ^ ) and is negligible due to condition (TJ|. 
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rather than r c , namely, 

S B = vrTir (d - r ) . (13) 

Note that the boundary condition at r = r c written in terms of tq reads 

d , , 

= - 5 ■ ("J 

Equations f[TH]) and (|14p provide the solution for the exponential factor 
in the probability of the merger of two strings for arbitrary ratio of the 
tensions T^ji^Tx). In particular at T2 = the problem is equivalent to that 
of spontaneous reconnection of two parallel strings [6]. From (TTUT) we learn 
that in this limit 

r = r c , (15) 

and the bounce configuration is described by a configuration discussed in [7] . 
In the present paper we will focus on the opposite limit in which the binding 
of the strings is parametrically small, i.e. 

2T x -T 2 <Ti. (16) 

It is convenient to introduce a dimensionless small parameter 5 for the bind- 
ing, 



<^/l-J§- (17) 

One can readily verify that in the limit 5^1 the gradient of the deviation of 
the string profile from a flat string is small, and the equations for the profile 
of the string world sheet in the bounce configuration can be linearized. This 
allows one to consider the fusion problem in a more physical strip geometry, 
i.e. with the strip boundary conditions. This can be done both for parallel 
strings and slightly nonparallel ones. 

2.2 Linearizing the problem in the weak binding limit 

In the linearized approximation (valid if |Vz| <C 1) the classical equation of 
motion for the string profile is 

Az = 0. (18) 
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If 8 <C 1 the above condition is met. It is not difficult to solve Eq. (1181) 
with the round boundary condition. Alternatively, one can expand the full 
nonlinear solution in the disk geometry. One finds in the leading logarithmic 
in R approximation at 5 <C 1 

r =r c 5 and r * - . (19) 

Then the action Sb can be readily derived from Eq. (fT5|) . namely, 

s °"* T *2i£m{ 1 + °hm]}- < 20 > 

Please, remember that the exponent determining the decay rate is e~ Ss . The 
condition ([2]) justifies the quasiclassical approximation. 

It is interesting to note that the bounce action is mainly determined by 
the tension T and the distance d, rather than by the binding parameter 5. 
The formula (I20p also tells us that the bounce action becomes small, and the 
semiclassical treatment becomes inapplicable, for exponentially long strings. 
However it is clear from the overall proportionality of the fusion rate to the 
string length that for such long strings the probability of fusion becomes of 
order one. It can also be readily verified that introducing of a small mass 
fi for the three-string junction, neglected throughout this paper, does not 
change the infrared dependence of the bounce action. Indeed, the /i-induced 
contribution to the action is 

ASb = a ""-»- iwm • (21) 



which has the same logarithmic behavior at large R as Sb in Eq. fl20|) . Thus, 
the condition under which the junction mass can be neglected does not de- 
pend on the string length and reduces to 

H < d <57i . (22) 

2.3 Strip boundary condition in linear approximation 

In the linear approximation in which the equation of motion for z(x, r) re- 
duces to the two-dimensional Laplace equation (fTgj) . the solution can be 
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constructed as a real (or imaginary) part of a holomorphic function of the 
complex variable 

w = x + ir . (23) 

Using this construction and the analogy with two-dimensional electrostatics 
(the so-called image charges method, see Appendix B), one can readily find 
the solution for the strip —L/2 < x < L/2 with the boundary conditions 
z(±L/2) = d/2. For the bounce centered at x = I and r = the solution 
has the form 



z(x, t) = r Re 



In 



sm 



■7r(w — I) 

2L 



COS 



2L 



+ 



d 



(24) 



with the constant r being determined by the condition of equilibrium of 
the bifurcation boundary, corresponding to z — 0, similarly to Eq. ffTO]) . 
Clearly, in the strip geometry the 0(2) symmetry is lost and the world sheet 
boundary for the string T 2 is no longer a disk. However for large L and for 
the bounce center not too close to the strip edge, the exact solution fl2~4l) can 
be approximated by a logarithmic one, 

w — I 



r Re 



In 



L 



+ const 



(25) 



It corresponds to an approximately circular bifurcation boundary with the 
radius r c related to the parameter ro as in Eq. ffl9|) . The applicability condi- 
tions for this approximation are that I is not parametrically close to L/2 and 
also L 3> r c . The bounce action 5*^ on such configuration in the logarithmic 
in L approximation coincides with that in Eq. (TSUI) with R being replaced by 
L, 

d 2 



S 



B 



2 \n(L5/d) 



1 + 



\n(L/d) 



2T X - T 2 



(26) 



3 Nonparallel strings 

When the number of the space-time dimensions is 3 + 1 (or more), it is 
possible to have nonintersecting and nonparallel strings. Then geometry of 
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the problem can be characterized by two parameters: the minimal distance 
d, and an angle a, which we will assume to be small (Fig. We chose 




Figure 5: Geometry for nonpar allel strings. 

the spatial axes in the following way, the z axis runs along the common 
perpendicular to the strings, and the origin of the coordinates is placed in 
the middle of the segment of this perpendicular connecting the strings. The 
x axis runs along the bisector of the angle between the projection of the 
strings on the {x, y} plane, so that the strings are parametrized as (a <C 1) 

7 /r. aX 

z = d/2, y = tan — , 

Oi X 

z = -d/2, y = -tan — . (27) 

We will refer to this configuration as "twisted strings." In choosing the in- 
frared regularization we aim at logarithmic accuracy of the bounce action, so 
that it is sufficient to consider the axially symmetric geometry of the world 
sheet with a large radius R, i.e. the round boundary conditions (Fig. [3]). 
Namely, the following constraints will be imposed: 

z(R) = ± d - (28) 
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for two strings under consideration, while the boundary conditions for y are 

OL 

y(x,r) r=R = ±-x . (29) 

In the linearized approximation the equations of motion for the orthogonal 
deviations of the string z(x, r) and y(x, r) are independent from one another. 
If a is small, "twisting" in the variable y can be considered as a perturbation 
over the solution z(x, r) for parallel strings (Sect. [2]). It is clear from the 
symmetry of the problem that y must identically vanish in the central part 
of the bounce, i.e. at r < r c one has y — 0. Then between the circles at 
r = r c and r = R the function y is harmonic and changes from y(r c ) = 
to the values prescribed by the boundary conditions (1291) . Invoking the well- 
known central harmonics for the two-dimensional Laplace operator, we find 
the explicit form of the profile y(x, r) for the upper and the lower strings, 

* = ± iT=k*( 1 -'?)- (30) 

Thus, the twist in the y direction results in an additional contribution to 
the (linearized) bounce action which takes the form 

S y S B = -^T ir 2 c . (31) 

Using the relations (fT9l) for r c one finds that the fusion probability for twisted 
strings T(a) is reduced in comparison with the parallel strings, 

7tt a 2 „ d 2 



T(a) = T(0) exp 



64 5 2 Tl In 2 (Rd/d) 



(32) 



Here T(0) is the merger probability for the parallel strings. It can be noted 
that, although the contribution to the exponential factor associated with the 
twist is of a higher order in 1/ ln(i?<i _1 ) than that in T(0), it has a nontrivial 
singular dependence on 5. The latter dependence implies that the string 
merger at weak binding takes place only if the angle between them is also 
small, i.e. at a < a max where 

\n(R5/d) 

"max ~ 1= , • (33) 

y 1] d 
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4 Fusion of parallel domain walls 



We will analyze the domain wall fusion at 2T\ — T 2 <C Xi, when the linearized 
approximation is applicable. In contradistinction with the string problem, in 
the domain wall problem we do not need any infrared regularization, since the 
solution of the three-dimensional Laplace equation falls off as 1/r rather than 
logarithmically. The coordinates that parametrize the Euclidean 3-brane are 
x, y, r, while the walls are given by the height functions z = f(x, y, r). The 
boundary conditions are set at infinite r = [x 2 = y 2 + r 2 ) 1 / 2 , 

z(r) — > ± - at r — > 00. (34) 

The solutions of the linearized equation Az = for the top and bottom 
walls are 

, , Ad 

zi{x,y,T) = + -, 

r 2 

z 2 (x, y,r) = I , (35) 

where the bounce is assumed to be centered at the origin. The two T\ branes 
meet at r = r c where 

r. = ^. (36, 

It is obvious that r c is the radius of the word volume of the composite wall 
(i.e. the world volume radius of the T 2 brane configuration) at the moment 
it leaves Euclidean and enters the Minkowski space (r = 0). 
The total Euclidean action is the sum of two contributions 

47T f°° z' 2 

S = (T 2 - 2Ti) — r z c + 2Ti J Anr 2 dr — 

4:71 

= _(2T 1 -T 2 ) — rl + nnr.d 2 , (37) 

where the first one comes from the composite wall in the middle, while the 
second from two tails of elementary walls. The action fl37j) is regularized: we 
subtracted the contribution of two parallel undistorted walls (Fig. Eh). In 
deriving this action we used Eq. (13T)|) . 
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Please, note that the signs are opposite. The first term is negative, and 
dominant at large r c . The second is positive and dominant at small r c . The 
bounce solution which is at the tip of hill can be obtained by extremizing 
Eq. (137|) with respect to r c , 



d / Ti 



2 V 27i - T 2 



^,3 / ^1 



S B = -Txd\ . (38) 

3 V 27i - T 2 v ; 

The probability of the wall fusion per unit time and unit area is proportional 

to e' SB . 

Now we can check that the linearization approximation is valid. The 
necessary condition is \z r \ -C 1 which is equivalent to A/r^ <C 1. Equations 
(1551) and flU imply 



A d /27i - 
r c V T x 



~ — ~- i . (39) 



This condition is met at weak binding, i.e. 

J i 

Note that it does not depend on the inter wall distance d. However, the 
distance must be much larger than the wall thickness. It must be large 
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enough to ensure Sb 3> 1. In particular, the choice T\ d 3 ^> 1 does the job. 
Another condition that was assumed in the consideration above is that the 
tension of the three- wall junction (closed circles on Fig. [6fr) is negligible, so 
that the junction contribution to the action can be ignored. 



5 Parallel walls at strong binding 

If Eq. (|40j) is not satisfied, the wall binding is strong, and the bounce action 
must be treated beyond the linear approximation, 

Ayr f°° / , N 

Sb = y (T 2 - 221) rl + 2T X / 4vr r 2 dr [VTT^ 1 - 1 J . (41) 

It is not difficult to see that there exists an "integral of motion" analogous 
to (El), 



r 2 z' 



rl 



(42) 



The solution for the wall bounce must be such that the left-hand side of 

o 



is r independent. We use the notation r 2 , for this constant. Then the classical 
equation of motion reduces to 



,,: ±^J=, («) 



dr ^r 4 - 



implying the following relation between the parameters vq, r c and d: 

2 f°° dr d 

Another relation between tq and r c is the equilibrium condition coinciding 
with Eq. (Q. Considering the expression fj42|) for the integral of motion in 
the wall case, one readily finds from this condition that 

/ T 2 
A „4 / 1 1 2 



16 



Invoking Eqs. 
ed to 

Sb = 



and 



(2T\ 



T 2 — - 3 
; 3 



the bounce action in Eq. (j4ip can be transform- 
r 3 + 8tt Ti 



dr 



47T 

T 



— 7i 

3 1 




(46) 



The second term in the first line is an elliptic integral which we will replace 
by its value mandated by the condition (jUj) . The expression in the square 
brackets in the second line reduces to a combination of elementary functions. 
Indeed, 



dr 



3r 4 — Tq 



3r" 



C — Tn 



(47) 



V r4 ~ r o 

This can be directly verified by differentiating both sides in Eq. (|47p over r, 



As a result, the bounce action takes the form 



S 



B 



An 

T 



T 2 r c rf - 2 Ti r c * /r c 4 - r 4 + T x r' Q d 



(48) 



After substituting in the latter expression the relation (1451) one readily finds 
that the first two terms in parentheses cancel, and one is left with a simple 
formula for the bounce action in terms of ro and d, 

S B = —T l rld. (49) 

In order to express the parameter r in terms of the separation distance d and 
the wall tensions T\ and T 2 one needs to solve the transcendental equation 
f )44|) which results in an elliptic function. In the limit of weak coupling 
one recovers the results of the previous section, while in the extreme strong 
coupling limit, T 2 — >• 0, one has r c = ro and we find 

r(i) 



^0 



2^r(§) 



d = 0.381... d , 
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so that the bounce action can be written as 

2 

T x d 3 = 0.609... T x d 3 . (50) 

Thus, an estimate Sb = const. ttTi d 3 works in the entire range of 
variation of T 2 in the problem of the domain wall fusion: from T 2 = 2T\ 
(where 5 — > 0) down to T 2 = 0. 



£(!) 
r(f) 



6 An instructive example 

In this section we will consider a particular example in which weakly bound 
domain walls naturally appear. Supersymmetric Yang-Mills theory supports 
[8] critical domain walls whose tensions are exactly known [8,9]. Namely, the 
tension of the fc-wall is given by the formula 

T fc = iV 2 A 3 sin^ (51) 

for SU(iV) gauge group. Here A is a dynamical scale parameter. The tension 
of the composite k walls is less than k times the elementary wall tension (i.e. 
the k = 1 wall). 

Hence, at large N we have 

7T 3 

2T 1 -T 2 = — A 3 , T^vrA^A 3 . (52) 

Invoking Eq. fl38l) we conclude that the probability of two parallel wall fusion 
(per unit time per unit area) is proportional to 



cxp 



-^N'(Ad) 3 . (53) 



7 Conclusions 



In this paper we have considered a generic problem of restructuring (fusion) 
of solitonic objects, due to binding energy. We paid particular attention to 
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formulating, and solving, the problem in a generic way, without any reference 
to particular underlying theories, or mechanisms responsible for the soliton 
existence and binding (i.e. no reference to microscopic physics). In this way 
our results can be applied in every instance in which our assumptions (sum- 
marized at the very end of Introduction) are met. Although we have focused 
on strings and walls in 3 + 1 dimensions, our results can be easily general- 
ized to higher- dimensional branes and/or higher space-time dimensions. For 
example, if we consider two p branes in (p + 2)-dimensional space-time, the 
fusion probability per unite volume and unit time is 

f / rp \(p-l)/2 ] 

T P ~ exp I -C p+1 2~* (p - If [-^brj Pi d p+1 ) \ , (54) 

where Gt is the volume coefficient, 

Ct = ■ (55) 

The problem for strings in 3 + 1 dimensions has a peculiar infrared behav- 
ior, and requires a regularization. For domain walls, or higher-dimensional 
branes, no infrared regularization is needed. 

The study and computation of the fluctuations around the bounce so- 
lution is left for future work This problem is essential in order to compute 
the coefficient C in Eq. ([3]). In some circumstances, it is just a numerical 
coefficient, which have little impact on the physical behavior. In other cases 
it can be of crucial importance. Assume we want to study the problem of 
infinite, parallel strings. We thus take the problem considered in Sects. 12.11 
and 12.21 and send the cutoff R — > oo. In this case r c and r both go to zero, 
and the bounce action formally vanishes. The center of the bounce is located 
at the center of the circle of radius R. Translation of the center becomes 
approximately a flat direction, as R — > oo. This is a usual divergence that 
is absorbed in the volume dependence, so that in fact we compute the decay 
probability per unit time and unit length. But there is also another (nearly) 
flat direction - the one due to a rescaling (x, r) — > X(x,t). This is peculiar 
to the case of strings. In the weak binding limit, the right-hand side of ([7j) 
becomes (nearly) scale invariant. This extra (approximate) zero mode must 
be properly treated in calculating the factor C. 
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Is worth to mention a particular example of binding energy between 
strings. It arises in type IIB string theory where (p, q) strings are bound 
states of p Fl-strings and q Dl-strings. Networks of (p, g)-strings have been 
studied recently, since they naturally arise at the end of some string theory 
inflation scenarios [11]. Tunneling effects may be relevant for the evolution 
of these strings networks. 

Finally, we would like to mention that work on the D-brane fusion was 
carried out (e.g. [10]) in string theory, in a setting specific to string theory. 
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Appendix A: Connection with field theory 

In the case of weak binding, and parallel branes, the problem can be recast 
in a simple field-theoretic formulation. Let us discuss it for the specific case 
of parallel strings. The "true" action is the Nambu-Goto one, with a tension 
that depends upon the distance between the two strings. For small perpen- 
dicular fluctuations z, i.e. at weak binding, we can use x to parametrize the 
space coordinate on the world sheets of the strings, and the following action 
ensues: 

S = -T(z) J dxdr y/l -d^zd^z 

= J dxdr T{z) ^-1 + 1^9^ + ...^ , (A.l) 

where z is the distance between two strings, and T(z) is the "combined 
tension" as a function of the distance. If we define a scalar field 0, 

<f> = zy/f(zj, (A.2) 
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assuming that the z dependence of T is adiabatically slow, we get the action 
in the form 

S = J dxdr QfyW- V(0)J , (A.3) 

with the canonically normalized kinetic term and 

V{4>) = T{z{<P)) . (A.4) 

The shape of the effective potential V(<f>) is very similar to T(z), although 
the change of variables from z to 4> changes the functional dependence in the 
vicinity of the minimum at z — 0, Fig. [7J This effective potentail is of a 
"false plateau" type. It is flat almost everywere, apart from a small domain 
near zero, where it drops off by 2T\ — T^. 

Details of the potential shape near the minimum depend on microscopic 
physics that is responsible for the binding energy. The results presented in 
this paper do depend on these details. What was crucial was the fact that 
the effective potential is essentially constant and flat everywhere, and drops 
to zero at a very short distance (the string/wall thickness) from the origin. 



5 = 2Ti - T 2 



thick d z 




Figure 7: Potential as a function of the distance z. 

The flat plateaux means that we have a classical moduli space of vacua, 
which, in turn, corresponds to the fact that there are no long-range forces be- 
tween the two strings (or walls). The tunneling is similar to the conventional 
false vacuum decay studied for metastable potentials [4]. A key technical 
difference is that a barrier is present only in the form of kinetic energy. 

We can add the following regulator to the potential: 

V e = ez 2 {z-d) 2 , e^O. (A.5) 
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The total potential T(z) + V t (z) mantains z = as a true vacuum. The 
plateau is lifted, and z = d is a metastable vacuum. In this formulation, 
the decay probability is determined by the conventional bounce of the type 
relevant to the false vacuum decay. At e — > oo the thin wall approximation 
is valid. On the other hand, at e — > (the case we are interested in) we find 
ourselves completely outside of the thin wall approximation, albeit, as we see 
in the bulk of the paper, a bounce-like solution exists and can be explicitly 
found. 

The bounce becomes exceedingly shallower as we decrease e. The number 
of space-time dimensions is crucial here. In (1 + l)-dimensional theory, the 
solution asymptotically is logarithmic, and the boundary condition z\— z% = d 
can only be imposed at a finite distance. In three or more dimensions the 
fall-off is power-like, and the boundary condition can be imposed at infinity. 

Appendix B: Holomorphic potentials 

A good way to directly derive Eq. §M§ in Sect. 12.31 is provided by analogy 
with two-dimensional electrostatics through the well-known in electrostatics 
image charges method. To this end we use the trick of putting auxiliary 
"mirror" charges as shown in Fig. [HJ We extend the strip to the entire 
complex w = x+ir plane, put the original unit charge at w — I, and then add 
a series of '-1' mirror charges at = (2/c+l) L—l, and a series of '+1' mirror 
charges at w£ = 2kL+l, with k being any integer. (In fact k = in the latter 
case corresponds to the original charge whose potential is being calculated 
with the strip boundary conditions.) In this way we certainly satisfy, due 
to the symmetry of the system of mirror charges, the boundary conditions 
that z must vanish at Kew = ±L/2. We then can find the solution for the 
harmonic function z(x, r) in terms of the real part of a holomorphic potential 
$(«;) produced by the constructed system of charges: z = C\ Re$(u>) + C 2 
with C\ and C2 being arbitrary constants. The holomorphic potential from 
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2L 
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strip 



Fi gure 8: The system of mirror 'charges' relevant to deriving the solution in Sect. 12.31 
with strip boundary conditions. The shaded strip is the physical world sheet for the strings 
and all the sidebands are its 'images'. The filled circles stand for positive charges and the 
open circles are for the negative 'images'. 

the considered system of charges is given by 



Using the Euler's formula representing the sin function as a product, one can 
explicitly find the products in Eq. flB.ll) up to inessential (although infinite) 
multiplicative constant: 







2L 




(B.2) 
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From here the harmonic function (|24p ensues. It is certainly easy to check, 
a posteriori, that the function in Eg. (j24p satisfies the necessary boundary 
conditions in the discussed problem. 
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